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-- Double. Mesa Edited by Sandman on May 12, 1978 2:14 PM 

DIRECTORY 

DoubleDefs: FROM "doubledefs" USING [DAdd, DCompare, DSub], 

InlineDefs: FROM "inlinedefs" USING [ 

BITSHIFT, LongCARDINAL, LongDiv, LongDivMod, LongMuU], 

StringDefs: FROM "stringdefs" USING [AppendChar]; 

DEFINITIONS FROM InlineDefs, DoubleDefs; 

Double: PROGRAM IMPORTS StringDefs EXPORTS DoubleDefs ■ 
PUBLIC BEGIN 

DDivide: PROCEDURE [num, den: LongCARDINAL] 
RETURNS [quotient, remainder: LongCARDINAL] - 
BEGIN 

qq: CARDINAL; 
count: [0. .31); 
ITemp: LongCARDINAL; 
IF den.highbits = THEN 
BEGIN 
[quotient. highbits, qq] <- 

LongDivMod[[lowbits :num.highbits,highbits:0],den. lowbits]; 
[quotient . lowbits , remainder . lowbits] *- 

LongDivMod[[ lowbits: num. lowbits, highbits :qq], den. lowbits]; 
remainder .highbits <- 0; 
END 
ELSE 
BEGIN 
count ^ 0; 

quotient . highbits ♦- 0; 
ITemp <- den; 

WHILE ITemp. highbits # DO -- normalize 
ITemp. lowbits <- 

BITSHIFT[lTemp. lowbits, -1] + BITSHIFT[lTemp. highbits, 16]; 
ITemp, highbits ^ BITSHIFT[1 Temp. highbits, -1] ; 
count <- count + 1; 
ENDLOOP; 
qq *- LongDiv[num, ITemp. lowbits]; -- trial quotient 
qq <- BITSHIFT[qq, -count]; 

ITemp <- LongMult[den. lowbits, qq]; -- multiply by trial quotient 
ITemp. highbits ♦- 1 Temp .highbits + den.highbits*qq; 
UNTIL DCompare[lTemp, num] ^ greater DO 

-- decrease quotient until product is small enough 
ITemp <- DSub[lTemp,den]; 
qq 4- qq - 1; 
ENDLOOP; 
quotient, lowbits ♦- qq; 
remainder «- DSub[num, ITemp]; 
END; 
RETURN 
END; 

DMultiply: PROCEDURE [a,b: LongCARDINAL] 
RETURNS [product: LongCARDINAL] « 
BEGIN 

product *- LongMult[a. lowbits , b. lowbits]; 
product. highbits ^ 

product, highbits + a. lowbits*b. highbits + a. highbits*b. lowbits; 
RETURN 
END; 

DNeg: PROCEDURE [a: LongCARDINAL] RETURNS [LongCARDINAL] « 
BEGIN 

RETURN[DSub[[0,0],a]]; 
END; 

DInc: PROCEDURE [a: LongCARDINAL] RETURNS [LongCARDINAL] « 
BEGIN 

RETURN[DAdd[[l,0].a]] 
END; 

AppendDouble: PROCEDURE [s: STRING, a: LongCARDINAL] - 
BEGIN OPEN StringDefs; 

longZero: LongCARDINAL « [highbits:0, lowbits:0]; 
longlO: LongCARDINAL « [highbits:0, lowbits:10]; 
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xn: PROCEDURE « 
BEGIN 

charZero: CARDINAL « LOOPHOLE[ • 0]; 
r: LongCARDINAL; 
IF a ^ longZero THEN 
BEGIN 

[a, r] <- DDivide[a, longlO]; 
xn[]; 

AppendChar[s, LOOPHOLE[r .lowbits+charZero, CHARACTER]]; 
END; 
END; 
IF a « longZero THEN AppenclChar[s , '0] ELSE xn[]; 
RETURN 

END; --AppendDouble 
END. .. 



